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Abstract. When are quantum filters asymptotically independent of the initial state? 
We show that this is the case for absolutely continuous initial states when the quantum 
stochastic model satisfies an observability condition. When the initial system is finite 
dimensional, this condition can be verified explicitly in terms of a rank condition on 
the coefficients of the associated quantum stochastic differential equation. 



1. Introduction 

It is almost a tautology that laboratory measurements give rise to classical stochastic 
processes. For example, in quantum optics one usually detects, using a configuration of 
photodctcctors, the light of a laser which is scattered off a cloud of atoms, and the resulting 
photocurrcnt is a classical stochastic process [21 0]. It is subsequently of interest to infer 
as well as possible the state of the atoms from the observed photocurrent, which is the 
purpose of quantum filtering theory. This theory has been extensively investigated both in 
the mathematical literature|3] (see Ref. [3] for a recent review) and in the physics literature, 
where it is known under the name of quantum trajectory theory or the theory of stochastic 
master equations [11]. 

In order to implement the quantum filter, however, the underlying quantum model is 
presumed to be known. It is not evident, a priori, that good estimates will be obtained in 
the presence of modelling errors which are inevitable in practice. Questions of robustness 
to modelling errors are particularly subtle on a long time interval, and have received much 
attention in the classical nonlinear filtering literature, see, e.g., Ref. [8] and the references 
therein. In particular, asymptotic stability of the filter — the independence of the filter, after 
a long time interval, of the initial estimate of the system — has been shown to hold in a 
wide range of classical nonlinear filtering models, and is the starting point for more general 
robustness questions. The problem of asymptotic stability is related to the consistency of 
Bayes estimates and is of significant practical interest as it ensures optimal performance 
of the filter, after an initial transient, even under misspecification of the initial condition. 
To date, however, no such result is known in quantum filtering theory. 

The goal of this paper is to develop a criterion which ensures asymptotic stability of 
quantum filters. This observability condition for stability is a natural one: it is the require- 
ment that no two different initial states of the model give rise to an observation process 
with the same law. In the quantum optics example described above, this means that we 
must be able to determine precisely the initial state of the atoms if we have access to the 
full statistics of the photocurrent over the infinite time interval. If this is the case, then the 
filtered estimates of the atomic observables are insensitive to the initial state of the atoms 
after a long time interval, provided we restrict our attention to initial states that satisfy 
an absolute continuity condition. 

The basic method of proof is based on the classical counterpart of this result, which has 
recently been developed by the author [21], To extend this result to the quantum setting, 
it is most natural to work within an abstract quantum filtering setting which is a little 
more general than the usual setting in the quantum filtering literature [21 HJ. We set up 
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the problem in section [2] in the context of C*-algebraic Markov process theory in the spirit 
of Accardi, Frigerio and Lewis pQ. The proof of the main result can be found in section 
[3] In section |4] we elaborate on the absolute continuity condition required by our main 
result, and provide a simpler sufficient condition. In the last section [5] we investigate a 
class of quantum filtering models, defined through the solution of a Hudson-Parthasarathy 
type quantum stochastic differential equation with a finite dimensional initial system, which 
have important applications, e.g., in quantum optics. In this setting one may find explicitly 
computable rank conditions for the model to be observable in terms of the coefficients of 
the quantum stochastic differential equation and the observation model. 

Finally, let us note that the asymptotic stability of nonlinear filters is not only of interest 
by itself, but is also an important ingredient in the development of error bounds for filters 
under more general modelling errors or for approximate filters (see, e.g., Refs. [3|6] in the 
classical setting). In that case, however, it is typically necessary to obtain more quantitative 
bounds on the rate of stability. Let us also mention that observability, though sufficient, 
is not a necessary condition for stability. One could conjecture that a natural counterpart 
of the detectability condition in Ref. [21] is necessary and sufficient for the stability of 
quantum filters in the finite dimensional setting of section [51 as it is in the classical case. 

2. The quantum filtering model 

We will consider quantum filtering theory in the abstract setting of Feller-type quan- 
tum Markov processes in the spirit of Accardi, Frigerio and Lewis [Tl 110], One of the most 
important examples in practice is the quantum stochastic flow generated by a quantum sto- 
chastic differential equation with a finite-dimensional initial system; this particular setting 
will be investigated in detail in section [5j 

In this section, we introduce the quantum filtering model and fix the notation for the 
rest of the paper. Let us begin by defining the basic elements of the model. 

• A, the initial system, is a unital C* -algebra with state space S C A*; 

• {-Ft,* > 0} is a one parameter semigroup of contractive and completely positive 
linear maps from A to itself, with P a [X] = X V X G A and P t [I] = I Vi > 0; 

• 971, the universal algebra, is a Von Neumann algebra; 

• {97t t ] : t > 0} is a filtration of subalgebras of 97t such that (Ut>o^*])" = ^ an< ^ 
9Jt ] — A** (i.e., 97t ] is the enveloping algebra of A); 

• {$p : p E S} is a family of normal states on 971 such that the conditional expecta- 
tions $ p ( ■ |97T t ]) : 97t — ► 971 t ] exist for every t > and p G S. 

Remark 2.1. The requirement that A be unital is not overly restrictive; if A is not unital, we 
may always enlarge A by adjoining the identity without essentially changing the structure 
of the theory. When A is commutative, this corresponds to the one-point compactification 
of the spectrum (Ex. VII.8.5 in Ref. [7]). 

Before proceeding, we recall for the reader's convenience the definition of the conditional 
expectation in a Von Neumann algebra (see, e.g., Ref. [T5|). 

Definition 2.2 (Conditional expectation). Let 21, 2lo be Von Neumann algebras, 21q C 21 
and let $ be a normal state on 21. Suppose there exists a linear map <&( • |2lo) : 21 — > 2lo 
which satisfies the following properties: 

. $(/|2l ) = I; 

• $(A*A|2l ) > for all X G 21; 

• $(A*|2l ) = $(X|2l )* for all X G 21; 

• $(ArZ|2lo) = X<£>(Y\%o)Z for all Y G 21 and X, Z G 2l ; 

• $($(X|2lo)) = ${X) for all X G 21. 

Then <E>( • |2lo) is a conditional expectation from 21 onto 21q with respect to $. 
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It is not difficult to prove that any two maps P, Q : 21 — > 2lo which satisfy this definition 
are ^-indistinguishable, i.e., $(|P(X) - Q{X)\ 2 ) = (see, e.g., Thm. 3.16 in Ref. g]). 
Thus the conditional expectation, if it exists, is essentially unique. Existence, on the other 
hand, is not guaranteed in the noncommutative setting. 

We now return to our filtering setup. We will presume that there is a family {j t : t > 0} 
of * -isomorphisms j t : A — > 9Jt t ] such that the Markov property holds: 

$ p (j t+s (X)\Wl s] ) =j s (P t [X}) Vi,s >0, X€ A, p€S. 

Moreover, we presume that jo(A)" = SUl ] an d that 

$„{j {X)) = p(X) VXeA, P eS, 

i.e., the state p £ S can be interpreted as the initial state of the quantum Markov process 
jt ■ The latter plays the role of the signal process in classical filtering theory. 

Remark 2.3. In order that $> p (jo(X)) = p(X) for all p, it is necessary that every state p E S 
extends to a normal state on £EJl i . This forces us to work with the universal representation 
M 0] ~ A** as required above, see Thm. 1.17.2 in Ref. [18] . 

In addition, we must introduce the observations. To this end, we introduce the n- 
dimensional observation process {Y t k : t > 0, k — l,...,n}, where Y t fc is a self-adjoint 
operator affiliated to 2Jt t j and Yq = 0. Define 2Jt] to be the Von Neumann algebra generated 
by {Y k : < s < t, k = 1, . . . , n}. We presume that 

2) t] is commutative, j t (X) € 2)^ Vt > 0, X e A. 

The first condition is known as the self-nondemolition property, and ensures that the 
process {Y t } can be represented as a classical stochastic process (as is befitting of laboratory 
observations). The second condition is the nondemolition property, and ensures that the 
conditional expectations 7rf (X) := $ p (j t (X)|2J t ]) exist for every X G A and t > (see, e.g., 
Thm. 3.16 in Ref. [4]). The goal of the filtering problem is to compute these conditional 
expectations. This problem can be solved explicitly in specific models, as is known since 
the work of Bclavkin[3j; see Ref. [4] for an introduction and review. For the purpose of 
this paper, however, it will not be necessary to obtain explicit expressions for the filtered 
estimates tt^(X). 

Finally we introduce the following Feller-type assumption. We presume that for any 
choice of ii, . . . j tfc > and bounded continuous functions fi, ■ ■ ■ , fk : K — > M, 

$ P (A(Y tl ) ■ ■ ■ f k (Y tk )\m 0] ) = uzih, ..../;.../, f k )) 

for some Z(t\, . . . , t%, /i, . . . , fk) G A independent of p, and moreover 

* P (/i(y. +tl - y s ) ■ ■ ■ f k (Y s+tk - y s )\ti s] ) = j s (z(h, ..../,../, h)) 

for every s > 0. The latter assumption ensures, in a sense, that the observation process is 
time-homogeneous. An important example of a filtering model in which these constructions 
can be implemented is discussed in detail in section [5] 

The goal of the remainder of the paper is to study the dependence of the filter TTf(X) := 
$ p (jt(X)|2)t]) on the initial state p <S S as t — * oo. 

Definition 2.4 (Observability). Let 2J = (Ut>o ?)*])"■ The model is observable if there 
do not exist pi,p2 65 with p\ ^ p 2 and & Pl (Y) = & P2 (Y) for every Y & 2J. 

We will prove the following result. 

Theorem 2.5. If the model is observable, then 

$ pi (\^(X)-n?(X)\)^0 VXeA 

whenever the laws of the observations under <& Pl and $ P2 are absolutely continuous (i.e., 
if P is a projection in 2) and <j> P2 (P) = 0, then <j> pi (P) = 0). 
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We can obtain a sufficient condition for the absolute continuity of the observation laws, 
as required in theorem 12.51 in terms of the initial states. This is developed in section [4j 
In the finite-dimensional setting, discussed in section we will find explicitly computable 
conditions for the filtering model to be observable. 

3. Observability and filter stability 
The proof of the main result proceeds in two steps. First, we establish that 

<M|7Tf(X) -TTf (X)|)^0 

for X of the form Z(t±, . . . ,tk, fx, ■ ■ ■ , fk)- This holds without any further assumptions. 
Then, we show that the set of all such observables is total in A when the model is observable. 
A simple approximation argument then completes the proof. 

3.1. Stability of Z(t\, . . . , ifc, /i, . . . , /&). We begin by proving a simple lemma. This 
result is almost trivial — it is just the tower property of the conditional expectation — but 
one should verify that the conditional expectations do in fact exist. 

Lemma 3.1. For any p £ S, s > and t%, . . . , tk, fi, ■ ■ ■ , fk, 

nP(Z(t u . . . ,t fc ,/i, . . .,/*)) = - Y s) ■ --fkiYs+t, ~ 

up to $ -indistinguishability. 

Proof. First, note that by the nondemolition assumption 

%(fi(Y s+tl -Y.)-.. f k (Y s+tk - Y s )\Tl s] ) = j s (Z(h, ..../;,../, fk)) G %- 

Hence the conditional expectation with respect to 2) s i exists and 

7rP(Z(h, ...,t k ,fi,...,f k )) = $ p ($ p (fi (Y s+tl - Y„) ■ ■ ■ f k (Y s+tk - Y s )\m s] )W s] ). 

Moreover, the conditional expectation 

$ P (h(Y s+tl - Y„) ■ ■■ f k {Y s+th - Y a )\fQ a] ) 

exists as fi(Y s+tl — Y s ) ■ ■ ■ fk(Y s+tk — Y s ) £ 2)^, (this follows directly as 2} is commutative). 
Finally, note that observables of the form X fi(Y s+tl — Y s ) ■ ■ ■ fk(Y s+tk — Y s ) with X E 2) s i 
are weak* total in 2). Hence the maps $ p ( • |2) s] ) : 2} -> 2) s] and * p (* p ( • |9Jl s] )|2) s] ) : 2J -» 
2) s ] ar e both well defined. It remains to note that both these maps satisfy the definition 
of the conditional expectation. □ 

We can now prove the stability of Z(ti, . . . /i, . . . , fk)- By virtue of the previous 
lemma the setting is essentially classical (as all the objects involved live in the commutative 
algebra 2)), and we will exploit this fact explicitly in the proof. 

Proposition 3.2. Suppose that the law of the observations under $ P1 is absolutely con- 
tinuous with respect to the law of the observations under $ P2 . Then 

$ P1 (K 1 (Z(h, ....'/.../, f k )) - Tif (z(t l7 ....'/.../, A))|) ^> o 

for any ti, . . . , t k > and bounded continuous functions /i, . . . , /fc. 

Proof. We work exclusively on the commutative algebra 2). By the spectral theorem (Prop. 
1.18.1 in Ref. [IB]), there exists a measure space (il, J^", A) which admits a surjectivc *- 
isomorphism i : 2) — > L°°(f2, J^", A), and every state Q v induces a probability measure 
P v on f2 such that $ V (X) = ~E v (i(X)) for all X € 2). Moreover, there exists a classical 
stochastic process {y^ : t > 0, fc = 1, . . . , n} on such that 

t(/(y£ , . . . , Y t k /)) = f(y% , . . . , y\l ) V bounded measurable / : R e ^ R, 
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and it is straightforward to verify that 

L(%(f 1 (Y s+tl -Y s )---f k (Y s+tk -Y s )\Z) s] )) 

= E ¥J (/ 1 (y s+fl - y 8 ) ■ ■ ■ fk{y s +t k - y s )\%) 

where ^ = a{y r : < r < s}. Evidently it suffices to prove that 

E pi (|E pi -E p2 (&|9J)|)±^o 

whenever £ s = fi(y s+tl — y 8 ) • • • fk(y s +t k - Us) for all s > and pi,p% £ S which give rise 
to absolutely continuous observation laws. 

To proceed, note that by our absolute continuity assumption P pi \oy oa -c ~Pp 2 W OD - We 
can therefore apply the classical Bayes formula (Lem. 8.6.2 in Ref. [16]): 

E P2 (A|^)E P1 (6|^)=E P2 (A6|^) 
where A = dP pl \ &oa ,/dP P2 \& x = V t>0 &t)- Thus we find that 

E P2 (A|^)|E P1 (6|^)-E P2 (6|^)| - |E P2 ((A-E P2 (A|^))6|^)| P P2 -a.s. 
Taking the expectation with respect to P P2 , we obtain 

E P1 (|E P1 (6|^)-E P2 (C 4 |^)|)=E P2 (|E P2 ((A-E P2 (A|^))^|^)|). 
By Jensen's inequality 

E P2 (|E P2 ((A-E P2 (A|^))6|^)|)< XE P2 (|A-E P2 (A|^)|), 

where K = ||/i||oo • • • ||/fc||oo- But note that £j is measurable with respect to so 
by the martingale convergence theorem E P2 (A|^) — > A in L 1 (fi, P P2 ). Therefore 
E P1 (|E P1 (&|^) - E P2 (£ t |^)|) -> 0, and the proof is complete. □ 

Corollary 3.3. Denote by O C A the linear span of Z{t\, . . . ,tk, fi, ■ ■ ■ , /&) for all 

t\, . . . , tk, fi, ■ ■ ■ , fki and suppose that the law of the observations under $ P1 is absolutely 
continuous with respect to the law of the observations under <i> P2 . Then 

$ Pl (|7r t pl (Z)-7rf (Z)\) VZedO , 
where c\O denotes the (uniform) closure of O in A. 

Proof. Fix Z <G clO° and a sequence {Z n } C 0° such that \\Z n — Z\\ — >• as n — » oo. For 
every n < oo, we have ^ Pi (\ , k^ 1 (Z n ) — Tr^ 2 (Z n )\) — > as t — > oo; to see this, it suffices to use 
the linearity of the conditional expectation and the fact that the triangle inequality holds 
for | • | when we restrict our attention to a commutative algebra (i.e., | Ylii^i\ — Si 1-^*1 
provided that the Xi commute with each other and their adjoints). Reasoning in the same 
way, we find immediately that 

<M|7!f (Z)-tff (Z)|) 

< $ pi (\^(Z- Z n )\) +$ pi (|7rf(Z„) -*?(Z n )\) + <f> pi (K" 2 (Z n -Z)\). 
The first and the third term on the right are bounded above by \\Z n — Z\\. Hence 
hmsup$ pi (|7rf (Z) - Trf (Z)\) < 2 \\Z n - Z\\. 

t — J-OO 

The result follows by letting n — * oo. □ 
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3.2. Observability and approximation. From the previous corollary, we see that a 
sufficient condition for the stability of the filter is that cl O = A. We will show that this 
is the case if and only if the model is observable. In fact, we will prove a more general 
result, from which this statement follows. We begin with the following definition. 

Definition 3.4 (Observable space). For p\,p2 € S, we define the equivalence relation 
pi ^ pi whenever $ Pl (F) = $ P2 (Y) for every Y S 2). The Banach space 

O = {X £ A : pi(X) = Pi{X) for all p\,p% & S such that p\ ^ p 2 } 

is called the observable space of the model. 

The following result is key. 

Proposition 3.5. 0° is dense in O . 

Proof. Suppose that 0° is not dense in O. Then there must be an element X of O that 
is not in clO°. By the Hahn-Banach theorem, there exists an element tp E A* such 
that tp(Z) = for all Z e clO° and ip(X) ^ 0. Then either tp{X) + tp(X)* ^ 0, or 
i(ip(X) — tp(X)*) 7^ 0, so we may assume without loss of generality that ip is real-valued. 
In particular, we can write tp = tp\ — ip% where tp 1 ,tp 2 are nonnegative (e.g., Prop. 1.17.1 
in Ref. [IB])- But note that I £ O , so tpi(I) = <P2(I)- We can thus define pi,pi £ S by 
Pi = <Pi/pi(I) an d P2 = «P2/<P2(I), and we find that pi(X) ^ P2(X) and pi(Z) = p2(Z) 
for all Z <G clO . Now note that for any p £ S 

p(Z(h, ...,t k ,h,...,f k )) = <P P (h(Y tl )- ■ ■ f k {Y tk )). 

Hence we find that 

* P1 (A (Y tl )---fk(Y tk )) = $ P2 (/j (Y tl ) ■ • • f k (Y t J) 

for all ii, . . . , tk, fx , ■ ■ ■ , fk- As the set of observables of the form f\(Y tl ) ■ ■ ■ fk (Y tk ) is weak* 
total in 2), we conclude that $ P1 (Y) = $ P2 (Y) for all Y G 2). But then p 1 — p 2 , which 
implies pi(X) = p 2 (X), and we have a contradiction. □ 

We immediately find the following corollary. 

Corollary 3.6. Suppose that the law of the observations under 4> Pl is absolutely continuous 
with respect to the law of the observations under $ P2 . Then 

^> Pl (\n? 1 (X)-4 2 (X)\)^0 VIGO. 

Proof. Immediate from corollary 13.31 and cl0° — O. □ 

We may finally complete the proof of theorem 12.51 

Proof. (Theorem 12. 5p . The model is observable, by definition, if p\ ^ p2 implies pi = p%. 
Clearly this is the case if and only if O = A. The result follows directly. □ 

Remark 3.7. The proof of proposition 13.51 clarifies why it is important to work in the C*- 
algebraic setting, rather than starting off with an initial Von Neumann algebra. As the 
state space of a C* -algebra is dual to the algebra itself, we may employ the Hahn-Banach 
theorem as in the proof of proposition 13.51 to characterize the observable space. For a Von 
Neumann algebra, however, the space of normal states is prcdual to the algebra. To employ 
the technique used in the proof of proposition 13. 5( we would then have two options: we 
must either consider non-normal initial states, or prove density of O in O in the weak* 
topology on the initial Von Neumann algebra. The former is unphysical, while in the latter 
case corollary 13.31 can not be employed. It thus appears that the C* -algebraic setting is 
the natural setting in which our results can be developed. 
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4. Absolute continuity and randomization 



In our main result, theorem 12 . 51 we required that the initial state pi,p2 & S are such 
that the law of the observations under $ Pl is absolutely continuous with respect to the law 
of the observations under <j> P2 . One might expect that a sufficient condition would be that 
the initial states are themselves absolutely continuous in a suitable sense. The goal of this 
section is to develop this idea. 

Before we turn to the filtering model of section [2] let us consider the general setting 
where A is any unital C*-algcbra. Given a state ip on A, we denote by (tt v , H^,,^) the 
cyclic representation of A induced by (p. 

Let S C A* denote the state space of A. We endow A* with the weak* topology, and 
recall that this makes S a compact convex set. By a (finite) measure on S we mean a 
regular Borel measure on S or, cquivalcntly, an element of C(S)* (see p. 232 in Rcf. [19]). 
A probability measure is a nonnegative measure with unit mass. 

We now recall a basic construction in Choquet theory. Let /i be a probability measure 
on S. Then (Lem. IV. 6. 3 in Ref. [19]) there is a unique p e S such that 

F(p)= [ F{ V )p(dip) V F £ A**. 
Js 

The state p is called the barycenter of the probability measure p. The measure p can be 
thought of as a randomization of the state p\ indeed, we have replaced the state p by a 
random state, with law p, which averages to p: 

p(X)= [ <p(X)p(dtp) VXeA. 
Js 

The idea is now to seek randomizations which have desirable probabilistic properties. In 
particular, we will consider the following notion of absolute continuity. 

Definition 4.1 (Absolute continuity). The state p\ € S is absolutely continuous with 
respect to p 2 G S, denoted as p\ <C P2, if there exist probability measures pi, p-i on S such 
that p\ is the barycenter of p\, p 2 is the barycenter of p 2 , and pi <^ p 2 - 

We now show that this natural definition of absolute continuity of p\ with respect to 
p 2 is equivalent to the requirement that p\ is presque dominee (almost dominated) by pi 
in the sense of Dixmicr (Ch. I, §4, Ex. 8c in Ref. [9]). Radon-Nikodym type results in this 
setting have been investigated by Naudts[P3] and Gudder[12]. 

Proposition 4.2. Let px^pi G S. Then the following are equivalent: 

(1) Pl < P2 ; 

(2) For every sequence {X n } C A such that lim mi „ pi((X m — X n )*{X m — X n )) = 0, 
we have lim„ pi(X*X n ) = whenever lim„ p2(X*X n ) — 0. 

(3) There exists a positive self-adjoint operator T on H P2 , affiliated to Tr P2 (A)' , such 
that Pl (X) = (T£ P2 ,ir P2 (X)T£ P2 ) for all X e A. 

Proof. (1 2) As pi <C P 2, there are probability measures p, v on S with p <^v and 



Pi(X)= / <p(X)p(d<p), P2{X)= / <p(X)v(dip), VXeA 
Js Js 

Let {A n } be such that lim mi „ P \((X m — X n )*(X m — X„)) = 0, and define the random 
variables $„ : S — > [0, oo[ by & n (tp) = <f{X*X n ). We begin by establishing that {$„} is a 
Cauchy sequence in L 1 (5,/i). By Lem. 2.4 in Ref. [14] 

|*m(¥>) - *n(v)| < <P((X m - X n )*(X m - In)) 1 ' 2 \ip(X* n X m ) 1 ^ 2 + ifiiX*^) 1 / 
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Therefore, we find using the Cauchy-Schwarz inequality and (a + b) 2 < 2a 2 + 2b 2 



Js 



< 2pi(X?nX m + X* t X n ) pi((X m - X n )*(X m - X n )). 



Thus {$„} is Cauchy in L l (S, p) provided that pi(A*A n ) converges to a finite limit. 
To show this, define tp n e H P1 by ip n = Tr Pl (X n )t; pi . Then pi(X*X n ) = \\ijj n \\ 2 and 
Pi((A" m — X n )*(X m — ^n)) = \\ipm — 4>n\\ 2 - As the latter converges to zero, we see that 
{ipn} is a Cauchy sequence in H P1 and thus pi(X*X n ) has a finite limit. 

Now suppose that, in addition, lim„ p2(X*X n ) = 0. Then evidently $„ — > in L 1 ^, v), 
so that in particular $„ — > in ^-probability as well as in /i-probability (as p <C ^ ). But 
as {<&„} is a Cauchy sequence in L 1 (5,/i), it follows that <f>„ — > in L X (S, /u). Thus 
lim n pi(Af *X n ) = 0, which is what we set out to prove. 

(2 3) See Cor. 2 in Ref. [12]. 

(3 1) Denote by £ the commutative Von Neumann algebra generated by T (i.e., this 
is the smallest Von Neumann subalgebra of Q3(H P2 ) which contains the spectral projections 
of T). By §3.1 in Ref. [18] ; there is a unique probability measure v on S with barycenter 
P2 and surjective "-isomorphism r : £ — * L°°(S 1 v) so that 

(£ P2 ,ir P2 (X)V£ P2 ) = f T(V)(cp)cp(X)is(d<p) V V G <£, X £ A. 
Js 

Now define f n {x) = nx/(n+x) and set T„ = f n (T). Then T n G £ is a bounded, self-adjoint 
operator and, writing the spectral measure of T as £r(dA), we find that 

||T„U - n fl2 \\ 2 = [ |/„(A) - A| 2 (£ P2 ,£t(^A)£ P2 ) ^+ 

J [0,oo [ 

by dominated convergence (as |/„(A) - A| 2 < 2f n (X) 2 + 2A 2 < 4A 2 , and ||T£ P2 || 2 < oo by 
construction). Consequently, we obtain 

r(T„ 2 )(^(AX^) = (T^ P2 ,7T P2 (X)T^ P2 ) Pl (X). 



is 

But note that T(T 2 )(ip) is nonncgative, nondecreasing and 



rK)(^(#)™i, 

so by monotone convergence T(T 2 ) /* A with A € L X (S, v). Thus 
Px{X)= [ A((p)(p(X)v(dkp) VXeA 



by dominated convergence. We now define dp = Adv, and p has barycenter pi. □ 

We now return to the filtering setting of section [2j The following result establishes 
that absolute continuity of the initial states is indeed a sufficient condition for absolute 
continuity of the observation laws. Absolute continuity of the initial states is often easily 
verified, e.g., in the finite dimensional setting of section [5] 

Proposition 4.3. Suppose that p\ <C P2- Then the law of the observations under $ Pl is 
absolutely continuous with respect to the law of the observations under $ P2 . 



Proof. As in the proof of proposition 13. 2[ we begin by constructing a measure space 
(CI, J^, A), a surjective "-isomorphism t : 2) — ► L°°(fi, A), and a family of probabil- 
ity measures P v on such that $> V (X) = E v (t(X)) for all X G 2). 

As p\ <C P2 , there exist two probability measures p\ , p2 such that p\ is the barycenter 
of p\, p2 is the barycenter of p2, and p\ -C P2- We will utilize these measures to construct 
randomizations of the classical probability measures P Pl and P P2 . To this end, let 8$ be 
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the Borel cr-algebra on the state space S, and construct the enlarged probability space 
(f2, #, P 2 ) by setting Q = Q x S, # = & x 3§, and 



^ 2 (d<p) VAgj^. 



P 2 (A)= / / /A(^,^)P,(d W ) 

Moreover, we define dP\ = AdP 2 where A(ui,ip) = {dpi/ dp2){ l fi)- Then 
[F(w)Pi(du,dfp)= [ F{u)P Pi (du;), i = 1,2, 

for all bounded measurable functions i* 1 : f2 — > R. As by construction Pi <C P 2 , clearly 
the marginals satisfy P Pl <C P P2 also. The proof is complete. □ 

5. The finite dimensional case 

In this section, we consider a specific class of quantum filtering models which have 
important applications in quantum optics (see, e.g., Refs. [U [ID- 
Fix p, q G N and let H = C p ® l~, where V = r s (L 2 (R+) ® C q ) is the symmetric Fock 
space of multiplicity q. Thus p is the dimension of the initial system, while q is the noise 
dimension. We set A = M p (the *-algebra of p x p complex matrices), 971 = 55(H) = 
A® 2J(0- Moreover, recalling that the Fock space admits the natural tensor product 
structure r = l~ t ] ® l~[ t , we define the filtration of subalgebras 

Tl t] ={X®I:Xe A® »(r t] )}. 

Finally, we define the family of states & p = p® $y with $y(X) = (£, where £ is the 
vacuum vector in T. It is not difficult to verify that the conditional expectations $ p ( • \9Jlt]) 
exist in this setting; in fact, they are given explicitly as follows: 

H t] : C ® V t] -> C p ® T, H$1> := V <8 £[*, $ p (X|9Jt t] ) = I^Xi^ <g> I, 

where £[ t is the vacuum vector in l~[ t . See, e.g., Ref. [17] for further details. 

We now introduce, as usual, the canonical quantum noises Ai(t), A\(t), Kij(t), i,j = 
l,...,q on T (we will denote their ampliations to H by the same notation), and consider 
the Hudson-Parthasarathy quantum stochastic differential equation 

{go go 
- 5 t] )dk. l3 {t) +Y J UdA\{t) 

- J2 L t S ki dA t (t) L* k L k dt -iHdt\ U u U = I, 

i,k=l k=l ) 

where qo < q and Sij,Li,H £ A, H is self-adjoint, and J2ij &ij ® e i e j ^ s a unitary operator 
in M p ® M q [ti is the ith basis vector in C 9 ). Then this equation has a unique solution 
{U t : t > 0} such that U t is unitary for every t > (Thm. 27.8 in Ref. [17]). Moreover, if 
we define jt : A — ► SE7l t ] by jt(X) = U^(X ® I)Ut, then j 4 satisfies the quantum Markov 
property for the semigroup {P t : t > 0} generated by 

JSf [X] = lim Pf[X j~ X = X] + J |^XL fe - i^L fc AT - ^L£Z fc | , 

fc=i J 

see Cor. 27.10 in Ref. [T7]. As by construction & p (j (X)) = p(X) for any p G S, this model 
satisfies the requirements of section [5] 

It remains to introduce the observations. For sake of concreteness, we will consider 
in detail two common observation models: a one-dimensional homodyne detection model 
and a one-dimensional photon counting model. The generalization of these results to 
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other observation models and to higher dimensional observations is straightforward. Before 
proceeding, however, we prove the following simple lemma. 

Lemma 5.1. Let p\,p2 be states on M p which are defined by the density matrices Q\,Q2 
(i.e., pi{X) = Tx\QiX]). Then pi <C P2 if and only ifkevgi D ker£>2- 

Proof. We first prove that ker g\ D kcr Q2 implies p\ <C p2- Let us restrict g\ , Q2 to the 
subspace h = (kcrg2) ± - Note that 02 |h has full rank and hence is positive definite, so 
there is some e G ]0,1[ such that (v,Q2v) > e|H| 2 for all v G h. But the eigenvalues of 
Qi must be contained in [0, 1], so that (v, g\v) < \\v\\ 2 for any d 6 h. Thus we find that 
(v,Q2v) > e(v, Q\v) for all v G h, so evidently Q2 > eq\. But then g[ = {Q2 — egi)/{l — e) 
defines another state p[ on M p , and the measures p\ <C P2 where p\ = <5{ Pl } and P2 = 
£ ^{pi} + (1 — £ ^{p[} nave barycenters p\ and p2- 

It remains to prove that p\ <C P2 implies ker D kerg2- To this end, suppose there is 
a v G ker ^2 such that v ^ kergi. Then p2(vv*) = (v, Q2V) = while pi(vv*) = (v, Q\v) = 
|| (gi) 1 / 2 w|| 2 > 0, contradicting p\ P2 by proposition 14.21 □ 

Note that by proposition ^. 31 this lemma makes the absolute continuity condition on the 
observation laws easy to verify explicitly in the finite-dimensional setting. In particular, 
the condition always holds if Q2 has full rank. This is very convenient in practice: it means 
that if the model is observable, we can always obtain the correct filtered estimates as t — > 00 
even when the true initial state of the system is completely unknown by choosing an initial 
state for the filter of full rank. 

5.1. Homodyne detection. For homodyne detection, we consider the observations 

Y t = U t *{^¥j(A 1 (t)+A\(t)) + ^/l^{A q {t)+A\{t))}U u rj G ]0, 1], q Q < q; 

here r\ is the detection efficiency, and the qth quadrature plays the role of an independent 
corrupting noise (we allow qo = q if r\ — 1). The operators Y t are self-adjointQ and 
affiliated to 3Jl t ] . Before we can proceed, we must verify that the nondcmolition and self- 
nondcmolition properties hold, as well as the Feller property of section [2j 

Lemma 5.2. Denote by 3ti the Von Neumann algebra generated by 

{Z s := VviMs) +A\(s)) + ^T^(A q (s)+Al( S )) :s<t}. 

Then U^3t]U T = U*3 t] U t = 2)t] for every < t < T. 

Proof. Denote by U s .t (s < t) the solution of the Hudson-Parthasarathy equation for Ut 
with the initial condition U s = I. Then it is not difficult to verify that U s> tU r _ s = U r ^ for 
r < s < t, and that U Syt acts as the identity on T s ] (Thm. 2.3 in Ref. [5]). Thus U s j & (3 S ])', 
so that U^ t ]UT = UtUt T 3t]Ut,TU t = U^ t ]Ut- Finally, note that any spectral projection 
P of Y s (with s < t) can be written as U*QU S where Q is a spectral projection of Z s , so 
that P = UtQUt also. But the set of all such Q generate 3t] and the set of all such P 
generate 2) t] , so U?3t] U t =%)t\. □ 

Corollary 5.3. The self-nondemolition and nondemolition properties hold: 

2)t] is commutative, jt{X) £ !)[] Vf>0, XeA 

Proof. As 3t] is a commutative algebra and 2} t ] = U^ t ]Ut, evidently 2) t ] is commutative 
also. To prove the nondemolition property, fix X £ A and P G 2) t j. Then [j t (X),P] = 
[U;{X®I)U u U* t QU t ] = U t *[X®I,Q]U t = 0, as [3Jt 0] ,3t]] = 0. □ 



The field quadrature Ai (t) + Al (t) should be interpreted as the Stone generator of the appropriate 
Weyl operator 1 171 . This defines the correct domain for these operators on which they are self-adjoint. 
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Remark 5.4. By virtue of the nondemolition and self-nondemolition properties, the filtering 
problem is well-posed. In this setting, one can compute the filter explicitly as the solution 
of the following stochastic differential equation: 

dt4(X) = ^(^[X])dt + ^j{^(L* 1 X + XL 1 ) + L\) v p t (X)} dW t , 

where dW P t = dY t ~ jrji${Li +L*)dt and ir?(X) = p(X), see, e.g., sec. 5.2.4 in Ref. [2D]. 
However, we do not need this representation of the filter in this paper. 

We must still demonstrate the remaining requirement of section [2j 

Lemma 5.5. For any t\, . . . ,t k > and bounded continuous fi, ■ ■ ■ , fk : R — > R, 

(*ti) ' ' ' fk(Y tk )\TX 0] ) = jo(Z(t!, . . . ,i fc ,/i, . . . , /*)) 

for some Z(ti, . . . , t k , /i, ■ ■ • , fk) S A independent of p, and moreover 

$p(/i(F s+tl -¥„)■■■ f k (Y s+tk - Y s )\m s] ) = j a (Z(h, ...,t fc ,/i,..., f k )) 

for every s > 0. 

Proof. The first assertion is trivial in the current setting, as the conditional expectation 
$ p ( ■ |9Jto]) does not depend on p and any element of SDto] can be written as jo(X) = X® I 
for some X G A. To prove the second assertion, note that 

h(Y s+tl -¥,)■■■ fk(Y s+tk - Y s ) = U* T h{Z s+tl -Z s )--- f k (Z s+tk - Z S )U T , 

where T is chosen to be greater than max{s + ti : £ = 1, . . . , k}. But as U s , U* € £DT S ] , we 
find by the module property of the conditional expectation 

*p(/i(y.+ti - y.) ■ • ■ A(n+t» - n)|an s] ) 

= f7* $ p ([/ s * T / 1 (Z s+tl - Z s ) ■ ■ ■ f k (Z s+tk - Z s )U S!T \Wl s] ) U s . 
Now note that for any pair of exponential vectors e(/), e(g) £ T and v, w G C p 
(« (» e(/), Ul T h(Z s+tl - Z s ) ■ ■ ■ f k (Z s+tk - Z S )U S , T w ® e(. 9 )) 

= (v®e(O a f),UZ_ B fi(Z tl )---f k (Z th )UT-.w<ae(e a g)) (e(f s] ),e(g s] )) 
= (v®e{6 s f),h{Y tl )---f k {Y tk )w®e{e s g)) (e(f s] ),e(g s] )), 

where s f(t) — f(s + t) and / s i is the restriction of / to [0, s]. Hence 

(v ® e(/), ^(C^T/i^+t! - Z s ) ■ ■ ■ f k (Z s+tk - Z s )U s<T \m s] ) w ® e(g)) 
= (v® e(f s] ) ® f [a> Ul T h(Z s+tl - Z s ) 

■ ■ ■ fk{Z s +t k - Z s )U s ,t w ® e(ff a] ) ® C[ s > (e(/[ s ), e(ff [s )) 
= <«®e,/i(^ 1 )---/fc(i r tj«'®0 (e(/ s] ),e( 5s] )) (e(/ [s ),e( 5[s )) 
= (« ® e, /i(^) • • • /fc(FtJ «; ® (e(/), e(. 9 )) 
= (t; ® e(/), ® p (fi(Y tl ) ■ ■ ■ f k (Y tk )\m 0] ) w ® e( 5 )). 

The result now follows as the exponential vectors are total in V. □ 

We have now completed verifying that all the requirements of section [2] are met, and 
thus theorem [23] applies. The remainder of this section is devoted to the following problem: 
can one determine directly whether the model is observable on the basis of the coefficients 
Sij, Li, HI We will find that this is indeed the case, and we will give an explicit algorithm 
to test observability. Most of the work consists of the computation of the characteristic 
function of the finite-dimensional distributions of the observation process; we employ for 
this purpose a technique used by Barchielli[2j. 
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Lemma 5.6. For any = to < t\ <t2<---<tk, we define 

T tl tfc (Ai, . . . , Afe) = ^(e^i^^-^-O+iA?^-^,)}]^). 

TTien we can write 

T tl ,..., 4fc (Ai, . . . , A fc ) = e^^ 1 ^^)* 1 • • • e^+^v^H**-**-!)/, 

w/iere JT[X] = Z/fX + AXi. 

Proof. Let k : [0, oo[ — > R be locally bounded and measurable and define 

3 t («) = f/ t * exp ^ k(s) dZ s + ^J^ k(s) 2 ds^j U t . 
Using the quantum Ito rules, we find that 

90 

dHt(K) = iK(i)v^Ht( K ) £{j t (S lfc )<M fc (i) + jt(Su)d4(t)} 
fc=i 

+ Ht)Vvjt(Li + LI) Et(n) dt + i K (t)y/l - fj5t(«) (dA,(t) + dAt(t)). 
Similarly, we find that 

90 

<&P0 = it(JSf [Jf]) dt + ]T - dA«(t) 



i,fe=l 

Using the quantum Ito rules once more and retaining only the time integrals, 



£ {j t (S* ki [X,L k ])dA\(t) +j t ([Lt,X]S ki )dA i (t)} . 

,k=l 

lore and retaining only the time integrals, 

j t (X)~ t (K) = X + [ {j s (^[X]) + iK(s)y/^jj s (J(r[X])}E s (K) ds + martingales. 
Jo 

Thus evidently, if we define T t (K,X) = $p(j t (X)Z t (n)\9Jl 0] ), then 

j t r t ( K , x) = r t (K, se[x\ + iK(t),/rjje[X\). 

The result now follows directly by setting 

k{s) = Ai/[ 0itl [(s) + A 2 /[ tl)t3 [(s) + ■ ■ • + A*/[ tfc _ litfc [(s), 
then solving the equation for T 4 (k, X) with X = I. □ 
Proposition 5.7. The observable space O can be characterized as 

O = span{Jz? Cl jr dl Jz? C2 • •■ &° h Jtr dk I : k,c t ,d t > 0}. 

In particular, O is the smallest linear subspace of A that contains I and is invariant under 
the action of Jzf and J(f . The model is observable if and only if dim O = p 2 . 

Proof. First, we claim that O coincides with the linear span of T tlj tk (Ai, . . . , A^) for all 
t\, . . . , ife, Ai, . . . , Afc. To sec this, note that the characteristic function of the joint distribu- 
tion of Y tl , . . , , Y tk under the state $ p is precisely p(T tli ... jtfc (Ai, . . . , A&)) (up to a constant 
factor). As the finite dimensional distributions determine the law of the observations, we 
have pi ^ p2 if and only if 

Pi(Tt 1 ,...,t fc (Ai, . . . , A*;)) = p2{Y tl ,...,t k {Xi, . . . , Aft)) Vti, . . . ,tfc, Ai, . . . , Afc. 

Thus evidently every element of the linear span of T tl tk (Ai, . . . , A^) is in O. Conversely, 
suppose that X G O is not in the linear span of T tl; ... jt( . (Ai, . . . , Afc); then there must 
exist an Y € M p such that Tr[F*A] ^ 0, but Tr[Y*Z] = for all Z in the linear span of 
Tf*i,...,t fc (Ai, . . . , A fe ). But writing Y as a(g a - g b ) + if3(g c - Qd) with a,j3 £ R and Q a ,...,Qd 
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density matrices corresponding to states p a , . . . , pd, we find that either p a (X) ^= pb(X), or 
p c (X) pd(X), while nonetheless p a ^ pb and p c ^ pd- We thus have a contradiction, and 
the claim is established. 

We now claim that the linear span of T tl ,....t k (Ai, . . . , A&) coincides with the linear 
span of «5f Cl J^ dl ■ • • J£ Ck Jtf dk I . First, note that any element of the latter form can be 

obtained from T tl tk (Ai, . . . , A^) by taking derivatives with respect to t, and A^. This 

means, in particular, that any clement of the latter form is in the closure of the linear 
span of Tti,...,t h (Ai, ■ ■ ■ , Afc). But we are working in finite dimensions, so the linear span 
is already closed. It remains to show that any T tl: .... ffc (Ai, . . . , A&) is in the linear span 
of elements of the form «Sf Cl <T dl ■ ■■^ Ck Jf dk I. This is an immediate consequence of the 
Cayley-Hamilton theorem, and the claim is established. 

Finally, we must show that O is the smallest linear subspacc of A that contains / and 
is invariant under the action of and Jff. Note that O is clearly invariant under Jzf and 
Jtf and contains /, so the smallest linear subspace such that this holds is contained in O. 
Conversely, any element of O can be generated by applying «if and J(f to / finitely many 
times and taking finitely many linear combinations, and the smallest subspace must contain 
at least these elements. This establishes the claim. Note that the model is observable if 
and only if O = A, which is clearly equivalent to dimO = p 2 . The proof is complete. □ 

Using this characterization of O we can construct and explicit algorithm for verifying 
observability. To this end, we define the linear spaces Z n C A by 

Z = span{7}, Z n = span{Z„_i, 3?[Z n -i], J^[Z n _i]}, n > 1. 

Clearly every element of O will be in Z n for some finite n. Moreover, if Z n = Z n+ \ for 
some n = m, then it is true for all n > m, and in particular Z m — O. But this will always 
be the case for some finite n: after all, the linear spaces Z n grow with n, but dimZ„ can 
not exceed p 2 . Hence this construction is guaranteed to yield O in a finite number of steps. 
To implement the procedure as a computational algorithm, one could start with {/} in the 
first step, then apply the Gram-Schmidt procedure at every iteration n to obtain a basis 
for Z n . 

5.2. Photon counting. In the photon counting case, we consider the observations 

Y t - U;{r,A n (t) + (1 - t?) A qq (t) + Vv^~v) (Ai,(t) + A ql (t))}U u 

where r\ G ]0, 1] is again the detection efficiency and qo < q. Once again Y t is self-adjoint 
and affiliated to 9Jl t ], and we must verify the various properties of section [2J The proofs 
of these properties are identical, however, to the homodyne case, so there is no need to 
repeat them. We only collect here the required facts. 

Lemma 5.8. Denote by 9T t j the Von Neumann algebra generated by 

{N s := 77 A u (s) + (1 - ?y) A g? (s) + v 7 '?^ - v) (Ai,(s) + A ql {s)) :s<t}. 

Then U^K^JJt = UfyiqUt = 2)t] for every < t < T. In particular, the self-nondemolition 
and nondemolition properties hold: 

2)t] is commutative, jt(X) G 2)t] Vt > 0, X G A. 

Moreover, for any t\, . . . ,tk > and bounded continuous fi, ■ ■ ■ , fk : K — > K, 

$ P (fl(Y t J ■ ■ ■ JlSKo]) = MZ(h, .... /;,../, fk)) 

for some Z{t\, . . . , tk, /i, ■ ■ • , fk) G A independent of p, and moreover 

® p (fi(Y s+tl - Y s ) ■ ■ ■ f k {Y s+tk - Y s )\m s] ) = j s (Z(h, ...,t k ,fi,..., f k )) 

for every s > 0. 
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Proof. The proofs of these facts are identical to the proofs of lemma 15721 corollary [531 an d 
lemma [531 and are thus omitted here. □ 



Remark 5.9. Also in this setting one can compute the filter explicitly as the solution of a 
stochastic differential equation driven by the observations: 



dn%(X) =irP(£'{X})dt + 



-<{x) 



(dY t — rj tx^{L\L\) dt), 



Trf(LiLi) 

where tTq(X) = p(X). We will not need this representation of the filter in this paper. 

To proceed, we must adapt lemma [5751 to the current setting. 
Lemma 5.10. For any = to < t\ < ti < ■ ■ ■ < tfc, we define 

T tl ,..., tfe (Ai,...,A fe ) = $ p (e^= 1 ^^-^- 1 )>|OT ]). 

Then we can write 

T tl ,.„, tk (Ai, . . ., A fc ) = e (^ lAl -D^)*i . . . ^Me^-Dvjm-t^j^ 
where J[X\ = L\XL X . 

Proof. Let n : [0, oo[ — > R be locally bounded and measurable and define 

St(«) = U* exp (ij n(s)dN s ^J U t . 
Using the quantum Ito rules, we find that 

90 

i.j=i 

+ (1-7?) (e IK(t) - l)St(«) dA« fl (t) 

+ vKf^ 1 ' 1 - 1)S 4 (k) ^{^(5^)^(4) +i t (5i i )dA^(t)} 

j=i 

90 

+ » 7 (e < "W - ^^(^^{it^Lx)^^) + j t (L^ u )dA(t)} 

i=l 

90 

+ VvO~v)(e lK{t) - l)E t (K)J2{jt(Ll)dA q (t) + j t (L 1 )dA\(t)} 

j=i 

+ 7 7 (e^ t )-l)5 t («)it(^ii)& 
Using the quantum Ito rules once more and retaining only the time integrals, 

j t (X)E t ( K ) = X + f {j a (Sf[X]) + (e' iK(s) - l)vis(J[X])} S a («) ds + martingales. 
Jo 

Thus evidently, if we define T t (K,X) = $ p (j t (X)E t {K)\DJt 0] ), then 

| r t ( K ,x) = r t (K,JZ[x} + (e«« - \)n/[x\). 

The result now follows directly by setting 

k(s) = AiJ[o, tl [(s) + A 2 /[ tl , t2 [(s) + • • • + A fe J[ tfc _ ljtfc[ (s), 
then solving the equation for Tt(n,X) with X = I. □ 

The following result now follows precisely as before. 
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Proposition 5.11. The observable space O can be characterized as 

O = ^&n{^ Cl f dl ^ C2 ■■■3? Ck f dh I : k >Ci ,di > 0}. 

In particular, O is the smallest linear subspace of A that contains I and is invariant under 
the action of Jz? and ^ . The model is observable if and only if dim O = p 2 . 

Proof. The proof is identical to that of proposition 15.71 □ 

5.3. Some remarks. In this section, we have obtained precise characterizations of when 
a homodyne detection or photon counting model is observable (when the initial system is 
finite dimensional) . This yields a simple algorithm to test observability from which stability 
of the filter follows directly due to theorem 12.51 Even in the absence of observability, 
however, one can say something about the stability of certain observables using corollary 
13.61 The simplest such result is the following. 

Corollary 5.12. For the homodyne detection model (resp. photon counting model), the 
observable M = L\ + L\ (resp. L*Li) is always stable in the sense that 

<IV(|7rf (M)-Trf (M)|)-^0 

whenever the law of the observations under $ Pl is absolutely continuous with respect to the 
law of the observations under $ P2 . 

Proof. This is immediate from corollary 13. 61 and the fact that M = L\ + L\ = ,J€\T\ (resp. 
M = L\L X = J\T\) is clearly in O. □ 

In the physics literature, the observable M in this corollary is sometimes called the mea- 
surement observable. The fact that the measurement observable is always stable regardless 
of any other properties of the model was established for the homodyne detection case in 
§5.3.2 of Ref. [20] using a different method. 

We conclude this section with an example that highlights the importance of the absolute 
continuity of the observations in our results. 

Example 5.13. We consider the homodyne detection model, and let us choose qo = 1, 
Sn = I, H = 0, and Li = F/2 with F = diag{l, 2, . . . ,p}. By the previous corollary, the 
measurement observable M = F is stable in the sense that 

<f Pl (\^(F)-n^(F)\)^0 

whenever the observations are absolutely continuous as required by theorem l2.5l It is easily 
verified, however, that any state p with density matrix of the form g = diag{0, . . . , 0, 1, 0, . . . , 0} 
is a fixed point for the filtering equation in remark HT4l Hence $> Pl (l^ 1 (F) — 7rf 2 (F)\) -/-> 
when pi , p2 are two different states of this form. Evidently the absolute continuity re- 
quirement is essential. We refer to Ref. [21] for a discussion of the connection between the 
weakening of the absolute continuity requirement and the notion of controllability in the 
classical setting. 
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